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ABSTRACT 

Context. Observations of water lines are a sensitive probe of the geometry, dynamics and chemical structure of dense molecular 
gas. The launch of Herschel with on board HIFI and PACS allow to probe the behaviour of multiple water lines with unprecedented 
sensitivity and resolution. 

Aims. We investigate the diagnostic value of specific water transitions in high-mass star-forming regions. As a test case, we apply our 
models to the AFGL 2591 region. 

Methods. A multi-zone escape probability method is used in two dimensions to calculate the radiative transfer. Similarities and 
differences of constant and jump abundance models are displayed, as well as when an outflow is incorporated. 

Results. In general, for models with a constant water abundance, the ground state lines, i.e., lio-loi> lii-Ooo> and 2i2-loi, are predicted 
in absorption, all the others in emission. This behaviour changes for models with a water abundance jump profile in that the line 
profiles for jumps by a factor of ~ 10-100 are similar to the line shapes in the constant abundance models, whereas larger jumps lead 
to emission profiles. Asymmetric line profiles are found for models with a cavity outflow and depend on the inclination angle. Models 
with an outflow cavity are favoured to reproduce the SWAS observations of the lio-loi ground-state transition. PACS spectra will tell 
us about the geometry of these regions, both through the continuum and through the lines. 

Conclusions. It is found that the low-lying transitions of water are sensitive to outflow features, and represent the excitation conditions 
in the outer regions. High-lying transitions are more sensitive to the adopted density and temperature distribution which probe the inner 
excitation conditions. The Herschel mission will thus be very helpful to constrain the physical and chemical structure of high-mass 
star-forming regions such as AFGL 2591. 

Key words. ISM: molecules - radiative transfer - stars: individual: AFGL 2591 



1. Introduction 

' High-mass stars, more than their low-mass counterparts, play a 
critical role in the energetics and dynamical feedback into the 
interstellar medium (ISM). Despite this, the formation mecha- 
nisms of massive stars are still a puzzle, owing to their occur- 
rence in more distant crowded stellar clusters, shorter forma- 
tion timescales, and formation in regions of high visual extinc- 
tion. Despite observational obstructions, promising theoretical 
work has been done o n the format i on of high- mass stars (e.g. , 
Bonnelletal. 1998; K lessenll200ll iMcKee & Tan 2002, |2003t 
Klessen et al.ll2007 ). Whereas low- and intermediate-mass stars 
(M< IOM0) are formed through gravitational c ollapse and sub- 
sequent disk accretion (iPalla & Stahleij [19931) . radiation pres- 
sure on the dust in the infalling gas prohibits high-mass stars 
to sufficientl y accrete gas and thus limits the stellar mass at 
20 -40 Mo (LY orke & Krugel [1977|; IWolfire & Cassinellilll987t 
iBeech & Mitalas 1994) . However, observations show that con- 
siderably more massive stars exist. Two formation scenarios 
have been proposed to work out this contradiction: non-spherical 
accretion (Yorke & Sonnhalter 2002) and merging of lower mass 
stars jBonnell et al.ll998tlBonnell & Batel2002,) . In order to dis- 
tinguish between these two scenarios detailed observations, sen- 
sitive to a range of size scales and physical conditions, are re- 
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quired. For this, interferometers can provide the spatial distribu- 
tion of dust and gas at the smallest size scales, while single dish 
telescopes equipped with large format, wide bandwidth, detec- 
tors are capable of simultaneously imaging a large number of 
molecules (albeit at lower spatial resolution). 

Up to date, there is growing evidence that disks are present 
around massive young stars as reports of direct imaging of disks 
in regions of high-mass star formation are increasing in number 
(e.g., Beltran et al. 2004; Zhane et al. 2002; Sandell et al. 200l 
Ivan der Tak & MentenI 120051; Ivan der Tak et al.i 2006). In addi- 
tion to direct imaging of disk-like structures, outflows provide 
indirect evidence for disks around massive youn g stars. In par- 
ticular, collimated outflows are a sign of accretion^B euther et al.l 
2002). Outflows are well known phenomena accompanied with 
sites of low-mass star formation, transporting angular momen- 
tum away from the star. The situation in the formation of high- 
mass stars is less clear, since spatial resolution has been lack- 
ing to resolve the outflows. However, over the last years, out- 
flows from high-mass stars h ave received inc r easing attention 
(see r ecent revie ws _by, e.g., IShepherdl[200l 12005^ ICesaronil 
120051; lArce et al.l l200'7h . It is found that the higher-mass ob- 
jects appear to continue smoothly the coiTelation found in T 
Tauri stars between outflow and accretion signatures. These re- 
sults suggest that the disk-outflow connection found in low- 
mass pre-main-sequence stars extends to more massive objects, 
and favour the accretion scenario. Besides this. iKrumholz et aTl 
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( l2005h show that outflows around massive stars lead to a signif- 
icant anisotropy in the stellar radiation field, hence reducing the 
radiation pressure experienced by the gas in the infalling enve- 
lope. As a result, a larger fraction of the material in the envelope 
can be used in the formation process of a massive star as opposed 
to models without outflows. 

The dynamical and physical evolution during star forma- 
tion is accompanied by strong chemical evolution. The cold- 
core phase, dominated by relatively simple molecular species, 
is transformed into a hotbed of more complex molecules af- 
ter the protostellar object heats up the suiTounding gas. Some 
of these molecules, e.g., H2, CO, and H2O, play an important 
role in regulating the temperature, hence pressure, via the pro- 
cess of heating and cooling through line absorption and emis- 
sion. iDotv & Neufeldl ( 1 1997 1) constructed models for the ther- 
mal balance, chemistry, continuum and line radiative transfer 
in dense molecular cloud cores with an embedded protostar. 
They find that due to temperature and density gradients, mo- 
lecular abundances are not constant within dense molecular 
cloud cores, i.e., freeze-out of water onto dust grains for tem- 
peratures below the sublimation temperature of ~100K and 
a high abundance of water in regions where the temperature 
is larger than lOOK. A lso, this behaviour has been found by 
ICeccarelli et aTl (I2OOOI) who modeled the structure of the col- 
lapsing envelope around the low-mass protostar IRAS 16293- 
2422. Besides this, outflows not only affect the physical state 
of the interstellar gas, but are also accompanied with chemi- 
cal changes. Observational estimates of molecular abundances in 
outflow regions have suggested that particular species, e.g., SiO, 
CH3OH, HC N, H2CO, SO, and SO2 are significantly enhanced 
(e.g., Bachil ler & Perez Guti errez 1997). In particular for water, 
the elevated gas temperatures found in the outflows increase the 
rate coefficients for the neutral-neutral reactions, leading to an 
enhanced water production rate. Consequently, the abundance 
of water shows significant variations up to a few orders of mag- 
nitude from one region to another This large variation in abun- 
dance makes water a powerful diagnostic of the physical struc- 
ture of the region and of the fundamental chemical interactions 
between the gas and the grains. Observations, together with the- 
oretical modeling of the available rotational and ro- vibrational 
lines from species and their isotopes, are neccesary to shed 
light on the density, temperature, and chemical structure over 
a wide range of conditions relevant for regions of star formation. 
Water in particular, with its many rotational lines, is very well 
suited to study the physical conditions in regions of low-mass 
jvan Kempen et al.li2007i) and high-mass (this paper) star forma- 
tion. 

In this paper, we model the excitation of water in the high- 
mass star-forming region AFGL 2591. The aim of our research 
is to make predictions for the line strengths and profiles of 
various ortho- and para-H20 (o- and P-H2O) transitions to be 
observed with the Heterodyne Instrument for the Far Infrared 
(HIFI) and the Photodetector Array Camera and Spectrometer 
(PACS) on board of Herschel, and to study how the results de- 
pend on the adopted geometry and chemical structure. Resulting 
line strengths for specific transitions are then compared to exist- 
ing Infrared Space Observatory (ISO) and Submillimeter Wave 
Astronomy Satellite (SWAS) observations. 

2. The high-mass star-forming region AFGL 2591 

While most massive stars form in clusters, GL2591, located 
in the Cygnus X region, provides one of the rare cases of 
a massive star forming in relative isolation. This allows to 



study the temperature, density, and velocity structure of the 
circumstellar envelope without confusion from nearby objects. 
Large columns of dust and gas toward GL2591 block our 
view of the stellar photosphere, making it invisible at opti- 
cal wavelengths, but results in bright infrared emission. The 
infrare d source is associated with a weak radio continuum 
source (Campbell" 1984 iTrinidad et al.ll2003h and with a bipo- 
lar molecular outflow larger than 1' in extent (|Lada et al.| [ l984l : 
Mit chell etal.lll992t iHasegawa & Mitchelllll995h . Observations 
of the molecular cloud complexes in the Cygnus X region in- 
dicate th at most objects are l ocated at the same distance of 
1.7 kpc (Schn eider et al.l 1200 6'). However, in the following we 
adopt a distance to the Cygnus X region of 1 kpc, fixing the 
lumin osity at ~2x 10"* Lp, to allow comparison with previous 
work (Ivan der Tak et al.lll999 ^). The results presented here de- 
pend weakly on distance, making the adopted distance inconse- 
quential. 

AFGL 2591 has been observed in w a ter lin es over a range 
of excitation conditions. iHelmich et al.l (Il996h report the de- 
tection of more than 30 lines within the bending vibration of 
water at 6 jjm u sing the Short Wavelength Spectrometer (SWS, 
de Graauw et al. 1996) on board ISO i Kessler et al.' 1996|) 
van Dishoeck & Helmich (1996,) and lBoonman & van DishoecS 
(|2003[) derive, from observed H2O spectra around 6nm, water 
abundances of ~ 10"^ up to 10"^ in hot-core regions. fSnell et all 
(2000) observed with the Submillimeter Wave Astronomy 
Satellite (SWAS) the ground-state transition of 0-H2O, thereby 
detecting the H2O emission from the more extended, cold gas. 
They find that the H2O abundances estimated for the hot core gas 
is at least 100 times larger than in the gas probed by SWAS, the 
latter being 6 x 10"'° to 10"^. Boonman et al. (2003) observed 
with the Short and Long Wavelength Spectrometer (SWS and 
LWS) on board of ISO and with SWAS H2O lines toward deeply 
embedded massive protostars, and find that ice evaporation in the 
warmer envelope and freeze-out in the cold outer parts together 
with pure gas-phase chemistry reproduces the H2O observations. 
This conclusion is strengthened by ground -based observations of 
the H^'^O isotope (I van der Tak et alj|2006l) . 

Despite previous observations of the regi on, the futur e 
launch of ESA's Herschel Space Observatory (Pilbratt" ' 2005h . 
with on board the HIFI will provide key information on the phys- 
ical and chemical conditions in molecular clouds by means of 
observing many water lines with higher sensitivity. Moreover, its 
high angular resolution allows to separate regions. HIFI, a het- 
erodyne receiver, is designed to provide continuous frequency 
coverage from 480 to 1250 GHz (Band 1-5), while Band 6 will 
cover the 1410- 1910GHz frequency interval. Going up in fre- 
quency, the beam size decreases from 39" to 13", and together 
with a spectral resolution of 0.2 -0.4 km s ' this instrument will 
be a powerful tool to probe the kinematics, i.e., infall and/or out- 
flow, and chemical complexity of the AFGL 2591 region. In ad- 
dition, PACS will take spectral images over the spectral band 
from 57 to 210;um. The main advantages of this instrument are 
the 9" resolution together with a 5 x 5 pixel imaging capability. 
Although the PACS spectral resolution of A/AA ~ 1500 is insuf- 
ficient to probe the kinematics in regions such as AFGL 2591, 
mapping of the H2O and continuum emission will provide sup- 
plementary information over a wavelength range inaccessible to 
HIFI. 



3. IVIodel set-up 

The temperature and density structure of the gas and dust 
in the envelope surrounding the young stellar object (YSO) 
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AFGL 2591 has been determined by Ivan der Tak et al.1 (Il999l 
l2000 h. In short, the resulting profiles are based on a study of 
single-dish submillimeter spectroscopy, combined with molecu- 
lar line emission at (sub)-millimeter wavelengths. The temper- 
ature structure of th e dust was modele d with the dust radiative 
transfer program of lEgan et alj (Il988l) . solving for the thermal 
balance of the grains as a function of distance to the star The 
resulting temperatures follow a r profile in the outer parts of 
the envelope, i.e., distances greater than 2-3 x 10"^ AU from the 
star, whereas in the inner envelope, the temperature gradient is 
steeper than r " The densi ty structure is then obtained with the 
use of the RATRAN code ( Hogerheiide & van der 1^ 120001 to 
model the observed molecular lines, i.e., CS, HCO^, HCN, and 
H2CO, using the temperature distribution calculated by the dust 
code. It is found that the density follows a power law of the form 
n = n(){rlr())^", with a~l-2. Dust emission profiles are consid- 
ered as well. In the following, results are presented for a -I 
and implications of an adopted a = 2 profile are elaborated on in 
Section|5] 

The H2O line emission in this work is calculated by ap- 
plication of the numerical code of Poelman & Spaans (2005, 
I2OO6 ) in which the radiative transfer is solved by means of 
a multi-zone escape probability method. The code, 'yS3D', is 
tested extensively against benc hmark problerns pres ented in 
lvanZadelhoffetal.1 (|2002) and v an der Tak et al.l (|2005'). The ra- 
diative transfer in this work includes rotational transitions of 
the ground-vibrational state up to an upper energy level of 
~ lOOOK. Vibrational pumping has been neglected as it is nu- 
merically challenging to include a substantial amount of ground- 
vibrational rotational levels, as well as levels of the vibrational 
mode in two dimensions. Spectroscopic and colli sional input is 
taken from the molecular database as presented in lSchoier et al.l 
( l2005h . Collisions between H2O and He, scaled to H2 to account 
for the differe nt reduced mass of the collisional partner, are taken 
into account (iGreen et~al ] |l993b . The adopted gas and dust tem- 
perature in our models ex ceeds 20 K at any poin t. Therefore, col- 
lisional rates presented in lDubernet et al.l (l2006l) for collisions of 
H2O with H2 in the temperature regime 5 - 20 K are not taken 
into account. The influence of dust radiation on the excitation 
of the w ater molecule is taken into ac count using grain proper- 
ties from lOssenkopf & HenningI ( Il994i) . Model 5. We adopt the 
expression for the ortho-to-para ratio (OPR) of H2O, in thermal 
equilibrium, defined by 

(2/o + 1) 2(27 + l)exp(- ^°«^-^'n 
(2/p + 1) 2(27 + l)exp(-^Sd^^) 

where /□ and /p are the total nuclear spin, corresponding to 
whether the hydrogen nuclear spins are parallel (/o = l, TT) 
or anti-parallel (/p - 0, ti)- The sum in the numerator (de- 
nominator) ex t ends o ver all ortho (para) levels (J,K^,Kcy, see 
lMunimaetaII(ll987h . 

Throughout this paper, a two-dimensional set-up is chosen 
in which the models consist of 78 logarithmically spaced shells 
with an outer radius of ~ 50 000 AU. A constant Doppler param- 
eter (1/e width) of 2kms"' is assumed for all the models. Fig.[T] 
plots the adopted temperature and density structure as a function 
of radius for the models with a -I and 2. It is seen that the a -I 
models tend to have lower density and temperature towards the 
center of the cloud. The gas and dust temperature is assumed to 
be coupled, due to the high densities reached throughout most 
of the cloud. The calculated level populations at each position in 
the cloud are then used to compute velocity profiles of various 
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Fig. 1: The distribution of density (left) and temperature of gas 
and dust (right) as a function of radius in the cloud is plotted. 
Solid lines mean a = 1, dotted a -2. 

transitions of o- and P-H2O using the program SKY that calcu- 
lates the sky brightness distribution, part of the RATRAN code. 
SKY has been tested against other codes. For transitions with a 
rest frequency within the frequency coverage of HIFI, the result- 
ing sky brightness distribution is convolved with the appropriate 
beam (depending on the frequency of the transition). In addition, 
predictions are made for the spectrum to be observed with PACS. 
Line strength s and profiles are co mpared with existing ISO ob- 
servations (B oonman et al.ll2b03h and therefore convolved with 
the appropriate ISO-LWS beam, between 66" and 78" in diam- 
eter, depending on the frequency. 

Due to the complex structure and spectroscopic properties 
of the water molecule, inverted transitions (masers) are likely 
to occur Maser action, a process of intense radiation confined 
within small an gular sizes w ith beaming angles in the range of 
10 ' - lO^^rad (lElitzurlll992 ). necessitates higher angular reso- 
lution than achieved with our adopted multi-zone escape proba- 
bility method. Moreover, a large negative optical depth, r, mag- 
nifies the exp(-T) term in the escape probability which can ham- 
per convergence. Therefore, we have artificially prevented the 
escape probabilities from exceeding exp(l), even when the level 
populations imply that - r > 1 . A maser transition becomes satu- 
rated in the regime that BJexp(-T) ~ F, with BJ the rate for stim- 
ulated emission and F the effective pump rate. We have checked 
that this limit is not reached for this transition, or other maser 
lines. The masers produced by our models are thus unsaturated, 
although this conclusion depends on the details of the adopted 
velocity field. As a consequence, preventing the escape probabil- 
ities from exceeding exp(l) has no consequences on the calcu- 
lated line intensities since non-saturated masers have negligible 
effect on the remaining level populations. Nevertheless, because 
of this limitation no predictions are made for maser action, e.g., 
22 GHz maser transition. 

4. Results 

4.1. Constant water abundance 

As a first s tep, we assume the H2O abundance to be uniformly 
distributed. ISnell e t al. (2000) find, using SWAS observations of 
the lio-loi ground-state transition at 557 GHz, that the abun- 
dance of 0-H2O relative to H2 in giant molecular cloud cores 
varies between ~ 10"'" and 10"^. In particular for AFGL 2591, 
a value of 6x10"'' is found. In view of this large observed 
spread, a constant H2O abundance with respect to H2, ranging 
from 10"'" to 10"^ is chosen. Resulting line profiles are plotted 
in Fig. |2] integrated brightness temperatures are listed in Table 
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Fig. 2: Line profiles of ortho- and para-HaO transitions for mod- 
els with a-\. The water abundance, i.e., 0-H2O + P-H2O, is uni- 
formly distributed across the region and ranges from 10"'*' to 
10"^ with 10-1° {solid), lO-** {dotted), 10"** {dash- dotted), 10"^ 
{dashed). 



|411 Beside Hj''© lines, multiple H^^O and transitions will 
be observed with HIFI and PACS for high-mass YSOs. Hence, 



model results in the case X(H20) is 10 and 10 are represen- 
tative for Hj^O and Hj^O, respectively. 

In general, the ground-state lines, i.e., lio-loi, lii-Ooo and 
2i2-loi, appear in absorption. All other lines appear in emis- 
sion. The optical depths in the lio-loi, lii-Ooo and 2i2-loi tran- 
sitions are factors of > 10 higher than the optical depths in the 
other lines, causing self-absorbed line profiles. When increas- 
ing the abundance by factors of 10, the line shapes transform 
from single-peaked emission into double-peaked self-absorbed 
profiles for the 2o2-lii, 22i-2i2, and 3()3-2i2 transitions. The 
strongest transitions are found to be those arising from the warm 
inner region, e.g., 3i2-3o3, 32i-3i2, and 3o3-2i2. 



4.2. Wafer abundance jump profile 

In this section we calculate the excitation of water, assuming that 
the H2O distribution follows a step function, i.e., a low abun- 
dance in the oute r regions and a higher abundance in the inner 
zones. lPotv et al.l (2006) have constructed models for the chemi- 
cal evolution in a collapsing massive young stellar object (YSO). 
These models include a realistic evolution of the central source, 
as well as ad/desorption of ices from grain mantles as the grains 
fall in from the cool exterior into the warm interior. They find 
that this approach reproduces the step-function water distribu- 
tion inferred observationally (v an der Tak et al. |'2006). For this 
reason, models with a water abundance jump profile are com- 
puted. The position in the cloud where the jump occurs cor- 
responds to the position where the sublimation temperature of 
water is reached, assumed to be 100 K in our models. The wa- 
ter abundance outside the sublimation radius is kept constant at 
10"^, whereas for regions with temperatures exceeding lOOK, 
the water abundance is increased from 10"^ to lO"** for the dif- 
ferent models. Resulting line profiles are plotted in Fig. [3] inte- 
grated brightness temperatures are listed in Table l4.2] 

It is seen that the line profiles of the ground-state lines, i.e., 
lio-loi, 2i2-loi, and ln-Ooo, are self-absorbed, similar to the 
constant abundance models. Line strengths and profiles of the 
ground-state transitions of o- and P-H2O, i.e., lio-loi, 2i2-loi 
and lii-Ooo, for models with a jump of a factor of 10-100 
are equivalent to the line strengths and profiles in case of the 
model with a constant water abundance of 10"''. Larger jumps, 
however, do influence the line strengths and profiles for these 



Table 1 : Integrated intensities, JrMBAv [Kkm s '], for models with 
a constant HtO abundance. 



Transition Frequency 



Water abundance 





[GHz] 


lo-'" 


10-' 


lo-** 


10-^ 


ortho-HaO transitions 


lio-loi 


557.936 


-l.l(-2) 


-0.2 


-1.2 


-1.4 


3l2-3o3 


1097.35 


3.3(-2) 


0.4 


4.1 


21.8 


3i2-22l 


1153.12 


6.0(-3) 


7.2(-2) 


0.9 


12.0 


321-3l2 


1162.91 


4.1(-2) 


0.5 


5.1 


34.6 


221-2l2 


1661.01 


2.0(-3) 


0.1 


1.1 


0.6 


2i2-loi 


1669.90 


-0.8 


-6.7 


-17.7 


-24.0 


3o3-2l2 


1716.77 


6.4(-2) 


0.7 


4.6 


5.0 


330-321 


2196.34 


1.1 (-4) 


1.4(-3) 


1.6(-2) 


0.2 


4i4-3o3 


2640.47 


6.6(-4) 


7.1(-3) 


5.2(-2) 


0.1 


221-110 


2773.97 


7.9(-4) 


9.6(-3) 


-2.0(-3) 


-0.7 


para-H20 transitions 


2ii-2o2 


752.029 


2.1(-2) 


0.3 


3.0 


19.1 


2o2-lii 


987.924 


3.2(-2) 


0.4 


4.3 


20.7 


lii-Ooo 


1113.34 


-0.2 


-2.0 


-7.1 


-10.1 


220-2ii 


1228.80 


2.0(-2) 


0.2 


2.8 


18.7 



NOTE: a(b) means a X 10* 
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Table 2: Integrated intensities, JTmbAv [Kkrns^'], for models with 
a H2O abundance jump profile. The water abundance outside the sub- 
limation radius is kept constant at 10"'. 
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Fig. 3: Line profiles of H2O for the jump abundance models with 
a -I. The water abundances for positions where r< lOOK is 
10"^, whereas for r> lOOK the abundances range from 10 
(solid), 10"^ (dotted), 10"'' (dash dot dot), 10"^ (dashed), lO"'* 
(dash dot). 



transitions. Thus, the ground-state lines are partly influenced 
by the warm inner component, and partly by the cold outer 
layer. Higher-lying transitions differ from the constant abun- 



Transition 


Frequency 




Inner water abundance 






[GHz] 


lo-** 


10-' 


10-* 


10-5 


10-* 


ortho-H20 transitions 


lio-loi 


557.936 


-0.2 


-0.2 


0.1 


1.2 


3.5 


3i2-3o3 


1097.35 


1.9 


6.4 


9.6 


13.2 


17.0 


3i2-22l 


1153.12 


0.4 


4.5 


17.3 


19.8 


20.5 


321-3l2 


1162.91 


2.3 


9.9 


17.0 


21.2 


25.9 


22l-2i2 


1661.01 


0.7 


2.0 


4.6 


9.3 


18.7 


2i2-loi 


1669.90 


-6.6 


-6.3 


-4.5 


1.7 


15.3 


3o3-2l2 


1716.77 


3.3 


6.3 


10.8 


13.2 


21.1 


330-321 


2196.34 


1.3(-2) 


7.7(-2) 


0.2 


0.3 


0.4 


4l4-3o3 


2640.47 


3.9(-2) 


8.1(-2) 


0.2 


0.3 


0.5 


221-110 


2773.76 


2.3(-2) 


6.0(-3) 9.0(-3) 


2.7(-2) 


0.2 


para-H20 transitions 


2ii-2o2 


752.029 


0.8 


2.5 


5.5 


13.4 


21.3 


2o2-lii 


987.924 


1.4 


5.7 


7.8 


8.6 


12.5 


lii-Ooo 


1113.34 


-1.9 


-1.9 


-1.2 


1.5 


7.4 


220-2ii 


1228.80 


0.8 


3.5 


5.3 


6.5 


11.2 



NOTE: a(b) means a X 10* 
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Fig. 4: Line profiles of two ortho-H20 transitions. The solid line 
represents the line profile for a jump with an abundance of lO-** 
where T<100K and 10-'' where T>100K. The dotted and 
dashed line profiles represent the constant water abundance of 
10-' and 10-^, respectively. 



dance model in that larger line strengths are found even for mod- 
erate jumps. This is due to the enhanced inner warm water abun- 
dance, the region that favours the excitation of the high-lying 
transitions. 

The peculiar appearance of the line profiles can be explained 
as a superposition of two constant models, i.e., a small, warm 
and abundant component together with a cold, large and low 
abundant component. In other words, a peak or drop in the line 
profile corresponds to a peak or drop in one of the transitions, 
when the contributing lines are optically thin. However, the line 
profile when one or both of the transitions are optically thick is 
not a simple superposition. In this case, photons coming from 
line center are lost and only those photons that escape in the line 
wings will contribute to the final line profile. This can be seen in 
Fig.|4]where we plot, for two different transitions, the resulting 
jump profiles with an inner water abundance of 10-^, together 
with the constant water abundance profiles of 10-' and 10-^. It 
is seen that for the 2i2-loi transition, the contribution from the 
optically thick inner 10"^ part is limited to those photons es- 
caping in the line wings. However, the resulting jump profile is 
dominated by the optically thin outer 10-' component. On the 
contrary, the jump profile for the 3 12-221 transition is most heav- 
ily influenced by the inner component. 
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One can wonder whether HIFI's future high spectral res- 
olution observations can distinguish between a constant or a 
jump profile water distribution. Since high-lying transitions, 
e.g., 3i2-3o3, 32i-3i2, 22i-2i2, more than the ground-state lines, 
are influenced by the warm inner water abundance, these lines 
are favoured to examine the occurence of an abundance jump. 
HIFI's high sensitivity allows to probe features in the line pro- 
file of ~ 1 K in Band 6, down to ~ 0.05-0.1 K in Band 1. This 
5cr detection is reached with an integration time of 1 hour for 
R- 10"*. Despite the differences in line strengths, a jump of a 
factor of ~ 10-100 leads to line shapes similar to the profiles in 
the constant abundance models. Hence, we do not expect HIFI 
to be able to distinguish between the two scenarios for jumps 
up to a factor of ~ 100. However, jumps of a factor of > 1000 
lead to line profiles more flattened around line center, with peaks 
in the line wings. This feature is not seen in the constant mod- 
els for abundances up to 10"^ for the high-lying transitions, and 
thus can be used to determine the water abundance structure. 
Nonetheless, to distinguish between constant and jump water 
abundance models is not as straightforward. It is likely that line 
profiles chan ge under the influence of inhomogeneity and veloc- 
ity gradients dPoeknan & Spaaris ''2006). Hence, one can easily 
misinterpret the observed line profiles. Multi-line observations 
in the future will help to further disentangle the events that are 
responsible for the observed line profile. 

4.3. Cavity 

We next consider cylindrically symmetric models in which an 
outflow cavity is incorporated. Observations confirm a power- 
ful bipol ar molecula r outflow associa ted with AFGL 2591 (e.g., 
iBally &L ada 19 831 iLadaet al.l ll984. ~ 1.5pc in extent. In ad- 
dition, v an der Tak et al. 1(12006^ constructed a representation of 
the geometry of AFGL 2591, based on millimeter interferomet- 
ric observations. They conclude that they are observing a young 
protostar with a disk close to face-on (inclination between 26- 
38 degrees). Therefore, models are constructed to study the ef- 
fect on the line strengths and profiles of a low-density outflow 
with an opening half-angle 6, as function of inclination angle. 
For this, a model with an outflow cavity with half-opening an- 
gle 6- 15° is computed. Models with larger opening angles are 
also computed and discussed in Section|5] The water abundance 
outside the cone is set at a constant value of lO"**, and the den- 



Table 3: Integrated intensities, JrMfiAv [Kkm s '], for a model with 
an outflow cavity with half-opening angle 6= 15°. 



Transition 


Frequency 




Inclination angle 






[GHz] 


0° 


20° 


45° 


70° 


90° 


ortho-HnO transitions 


lio-loi 


557.936 


8.7 


6.8 


3.1 


1.1 


0.7 


3i2-3o3 


1097.35 


32.9 


17.4 


3.7 


1.8 


1.6 


3l2-221 


1153.12 


5.5 


2.7 


0.5 


0.3 


0.2 


321-3|2 


1162.91 


36.6 


17.8 


3.4 


1.7 


1.5 


22l-2i2 


1661.01 


17.5 


7.0 


0.8 


0.4 


0.3 


2i2-loi 


1669.90 


105.8 


33.9 


-5.0 


-7.2 


-6.8 


3o3-2l2 


1716.17 


67.2 


25.3 


3.6 


1.8 


1.5 


330-321 


2196.34 


0.3 


0.3 


0.3 


0.2 


0.2 


4l4-3o3 


2640.47 


1.6 


1.6 


1.4 


1.1 


1.0 
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2773.97 


2.7 
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Fig. 5: Line profiles of H2O for a model with a = 1, a cavity with 
half-opening angle 0= 15° and an outflow velocity field with 
Vout = 12kms"'. The abundance outside the cavity has a con- 
stant value of 10"^, whereas inside the cavity the abundance is 
10"^. The different line profiles are the result of the different 
inclination angles. Inclination of 0° means face-on, as the line 
of sight is directed along the cavity; inclination of 90° is edge- 
on. IncUnation = 0° (solid), 20° (dotted), 45° (dash dot dot), 70° 
(dashed), 90° (dash dot). 
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sity and temperature structure follow Fig.[T] The outflow is char- 
acterised as a low density region with n(H2) = 5 x 10"^ cm"^, in 
agreement with the overal l density found thr ough hydrodynami- 
cal modeling of outflows dShang et al.ll2006l) . and a coupled gas 
and dust temperature of 300 K. This high gas and dust temper- 
ature i s reasonable as in ou tflows the temperatures reach a few 
100 K (iGiannini et al.l200ll) . Due to these high temperatures, en- 
hanced neutral-neutral rates lead to elevated H2O abundances. 
For this reason, a constant water abundance of 10"^ is chosen 
throughout the outflow. The outcome of models with a water 
abundance of 10"^ in the outflow is discussed in Section |5] The 
outflow velocity field is modeled as a function of radius, i.e., 
V oc R, reaching a maximum velocity of + 12 km s"' at the outer 
radius. Resulting line profiles are plotted in Fig. |5] integrated 
line strengths are listed in Table 14. 3] 

One notices the diff'erence in width of the line profiles for 
models with an inclination angle of 0°, compared to 90°. This is 
the result of the fraction of the outflow velocity field lying in the 
line of sight, i.e., face-on models tend to have a larger line pro- 
file width as opposed to edge-on models. It is seen from Table 
I4.3l that a decrease in inclination angle leads to an increase in the 
line strengths as the cavity is more uncovered. The outflow is 
best seen in the ground-state transitions for models with a small 
inclination angle, i.e., close to face-on. This is expected since 
the low density cone favours the emission from low-lying tran- 
sitions, in that the critical densities of the water lines are of the 
order of 10** cm"^. From Fig. |5] we see that the peak in the blue 
wing is weaker than the peak in the red wing of the line profile 
for the 1 10-1 01 ground-state transition, in contrast with the line 
profiles of the constant abundance and jump profile models. Note 
that for an inclination angle of 0°, the optical depth through the 
center, i.e., along the outflow, is larger than 1 for ground-state 
transition. It is found that models with an inclination angle of 
90° produce line profiles similar to the profiles in the case of a 
model with a constant water abundance of lO"**, except for the 
transitions seen with ISO. However, line strengths differ by a 
factor of a few. 

Clearly, HlFl's sensitivity can put constraints on the inclina- 
tion angle of the outflow, in that for small inclination angles, the 
blue wing is found to be less strong than the red wing. 

4.4. Disk 

Even though the formation of high-mass stars is not completely 
understood, more and more observational evidence is found that 
the crea tion originates in a si milar manner as its low-mass coun- 
terpart jCesaroni et al.|[T999 l). Therefore, we examine the influ- 
ence of an embedded disk in the center of the AFGL 2591 region 
on the line profiles and strengths. We consider the scenario of 
Section l4n for which we adopt a constant water abundance of 
10"^ throughout the region. 

The size of a circumstellar disk can be quantified as 
ra = 350 AU (r/0.007) (M*/20Mo)i/^ where y is the ratio of ro- 
tational to gravitational e nergy. For high-mass star-formation, 
y ~ 0.007 jPirogov et all |2003i). When adopting a mass of 
the central star toward AFGL 2591 of 16 M©, as inferred by 
Ivan der Tak & Menten (2005), a disk with radius ~300AU is 
found. We adopt an overall disk density of 10^ cm""*, and a wa- 
ter abundance of 10"^. The temperature of the gas and dust is 
assumed to be lOOK. Note that the characterization of the den- 
sity and temperature structure of the disk is far from accurate. 
However, the intent of this section is to investigate the influ- 
ence of a 'simple' disk model, not to model a realistic disk. The 
excitation of rotational, as well as ro-vibrational, water transi- 



Table 4: Integrated intensities, J^mb A v [Kkms '], 
for two different models, when convolved with the 
PACS beam. 



Transition 


Wavelength 
[jum] 


Model I 


Model 11 


ortho-H20 transitions 


330-221 


66.48 


0.3 


13.2 


32|-2l2 


75.43 


-1.7 


27.0 


423-3l2 


78.79 


0.6 


11.7 


5o5-4l4 


99.56 


1.5 


4.3 


221-110 


108.1 


-1.9 


25.4 


4l4-3o3 


113.5 


2.7 


17.6 


432-423 


121.7 


1.8 


7.6 


423 -4l4 


132.4 


2.2 


7.9 


3o3-2l2 


174.6 


6.9 


36.5 


2i2-loi 


179.5 


-23.1 


39.2 


para-HiO transitions 


331-220 


67.13 


0.2 


8.3 


322 -2 11 


89.99 


0.4 


7.6 


220-1 11 


100.9 


-2.2 


14.8 


3l3-2o2 


138.5 


2.5 


14.4 



NOTE: In Model I a constant water abundance of 
10"** is assumed. Model II describes AFGL 2591 as 
a region with a cavity outflow, with = 12 km s"' , 
seen face-on, i.e., inclination angle = 0°. To convert 
to Wcm"^ multiply /FmbAv by x 10"". 



tions in a circumstellar disk when the chemistry and tempera- 
ture structure is calcula ted self-consistently is under construction 
(iMeiierink et al.ll2007l) . 

It is found that the line profiles and strengths are insignifi- 
cantly influenced under all possible inclination angles (less than 
1 %). Since HIFI will be able to spectrally resolve disk struc- 
tures, the emission from our adopted disk model with a diameter 
up to a few hundred AU in size is overwhelmed by that of the 
envelope within the HlFl beams. However, observations show 
that more massive, elongated disks with a diameter up to few 
thousands of AU are likely to exist (iBeltran et al.ll2005l) . Due to 
high densities and large columns of dust in circumstellar disks, 
the stellar radiation field is attenuated. As a result the temper- 
ature of the gas and dust at radial dist ances of a more than a 
few 10 of AU drops below 100 K (see iGlassgold et al.ll2007h . 
Consequently, water will be frozen out on dust grains, reducing 
significantly the gas phase water abundance. When adopting the 
same disk density as in previous model, a water abundance value 
of 10 and a disk radius of 1500 AU, the line profiles differ only 
a few % from the models without a disk. 

4.5. Predictions for PACS 

Besides HlFl detections, observations by PACS will enlarge sig- 
nificantly our knowledge on high-mass star-forming regions. 
PACS will caiTy out spectral scans from 57 to 210/im, thereby 
providing information on the dust continuum emission at wave- 
lengths currently poorly characterized. However, the PACS data 
will be spectrally unresolved (A/AA ~ 1500), in contrast to HIFI 
observations. 

In this section, we compute the spectral output for PACS 
for two different scenarios. The first scenario (Model I) is for a 
model with a constant water abundance of 10"^, the second sce- 
nario (Model II) is for the model as described in Section [431 in 
which the outflow cavity is seen under an inclination angle of 0°. 
Resulting spectra are plotted in Fig.|6] integrated line strengths 
are listed in Table |4.5] 
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Fig. 6: Predictions for the spectrum to be observed with PACS 
for two different scenarios. Integrated intensities are plotted on 
top of the continuum background. The top panel denotes the 
spectrum for a model with a constant water abundance of 10"**. 
Bottom panel plots the spectrum when a cavity with an outflow 
velocity field, with Vom - 12kms"', is seen under an inclination 
angle of 0°. The Y-axis is in Kkms To convert to Wcm"^, 
multiply by x 10"^^. 



It is seen from Fig. |6] that the two models lead to rather dif- 
ferent output spectra. Strong lines are 220-I11, 221-I10, 4i4-3o3, 
3o3-2i2, and 2i2-loi. Depending on the model, the 32i-2i2, 220- 
lii, 221-1 10, and2i2-loi transitions are found either in emission 
or in absorption. One notices the change in the behaviour of 
the dust continuum emission between the two models. Hence, 
PACS spectra will tell us about the geometry of these regions, 
both through the continuum and through the lines. Note that the 
signal-to-noise (S/N) ratios for the different lines vary between 
30 and 3000 for Model I and between 25 and 300 for Model II. 

It is no w interesting to compare the results from 
iKaufman & N eufeld (1996) in which far-infrared water emis- 
sion from shock waves is presented. They find that the strongest 
transitions are the 2i2-lio, 3o3-2i2, 221 -1 10, 4i4-3o3, 3i3-2o2, lio- 



Table 5: Observed ISO and SWAS line intensi- 
ties [Kkms"'] for AFGL 2591, see Boonman et 
al. 2003 

Transition and Wavelength [GHz] 
2773.97 2640.47 2196.34 557.936 

22I-I1O 4i4-3o3 330-321 1 iQ-lpi 

<-0.5 <0.7 < 1.2 0.37 ±0.04 



loi, and 322-2ii lines. We find that these transitions are also the 
strongest in the case of models with an outflow cavity. 

5. Discussion 

5.1. Comparison to SWAS and ISO observations 

Observations by SWAS of the ground-state tran sition of 0-H2O 
at 55 7 GHz reveal an asymmetric line profile dBoonman et al.l 
l2003h . However, radiative transfer modeling of Une profiles in 
the massive young stellar object AFGL 2591 show that con- 
stant and jump water profile distributions do not fit this SWAS 
observation of the ground-state transition of 0-H2O, in that a 
symmetric line profile is found. In contrast to these models, an 
asymmetric line profile is derived for the ground-state transi- 
tion for models with an outflow cavity and an inclination an- 
gle between 0° and 30° . This agrees with the observations by 
Ivan der Tak et al.l (l2006l) . who conclude that the outflow is seen 
under an inclination angle of 26°- 38°. The integrated intensity 
of +7-9 Kkms"' reduces, after convolution with a beam size 
comparable to the size of the SWAS beam, i.e., 3.3' x4.5', to 
a value of +0.5 Kkms"'. This value is in agreement with the 
observed 0.37 ±0.04 Kkms"' strength. Note that this observed 
value is a lower limit, and that the actual value may be increased 
by a factor o f ~ 2, due to the prese nce of foreground clouds as 
discussed in lBoonman et al.l (|2003|) . 

Table 15.11 Usts the observed integrated lin e intensities for 
some selected ISO transitions, published in iBoonman et alj 
(|2003|) . The observed line fluxes [Wcm"^/im"'] have been con- 
verted into integrated line intensities [K km s" ' ] , by adopting the 
appropriate ISO beam size, to simplify comparison with the re- 
sults presented in this chapter For the 33o-32i and 4i4-3o3 transi- 
tions, the predicted integrated intensities are below the observed 
strengths for all the models by a factor of ~4, except for the 
models with an outflow that overproduce the 4i4-3o3 strength 
by a factor of ~2. The 221-I10 line is observed in absorption, 
whereas our models favour to produce this transition in emission. 
The model with a constant water abundance of 10"^ is favoured 
to fit this transition. 

Mode ling efforts of th i s pape r are similar to the results pre- 
sented in iBoonman et al.l (l2003h in that (i) the 4i4-3o3, 33o-32i 
and 3()3-2i2 transitions are always predicted to be in emission for 
models with a -I. (it) The 221-1 10 and 2i2-loi lines are found in 
emission for models with an outer abundance of 10"''. (Hi) Even 
though we tabulate in Section|4T|the resulting integrated intensi- 
ties up to a constant water abundance of 10"^, the trend is visible 
that models with a larger abundance would result in a deeper ab- 
sorption of the 221-1 10 transition than observed. Hence, models 
without freeze-out are not favoured to model this transition, (iv) 
IBoonman et al.l (l2003l) conclude that models with an inner water 
abundance of 2 x lO"'* and an outer water abundance in between 
10"" and 10"^ are favoured the match the observations for all 
the transitions. Our jump model with an inner water abundance 
of 10"'* produces integrated intensities that match closely these 
prefeiTed models. 
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Table 6: Integrated intensities, JTmbAv [Kkms"'], for models with 
a constant H2O abundance. The adopted water abundance is denoted 
on top of each column. The density and temperature distribution is 
described by a cf = 2 profile coefficient. 



Transition 


Frequency 10 


io-« 


io-« 


10-' 




[GHz] 










ortho-H20 transitions 


lio-loi 


557.936 


-6.0(-2) 


-0.6 


-1.6 


-2.1 


3l2-3o3 


1097.35 


5.6(-2) 


0.6 


3.5 


9.9 


3l2-221 


1153.12 


1.8(-2) 


0.2 


2.4 


15.4 


321-3l2 


1162.91 


0.1 


1.0 


7.0 


23.7 


22l-2i2 


1661.01 


-0.2 


-1.6 


-8.2 


-18.3 


2i2-loi 


1669.90 


-2.3 


-14.5 


-31.9 


-42.4 


3o3-2l2 


1716.77 


-0.1 


-1.1 


-6.4 


-16.8 


330-321 


2196.34 


-l.l(-3) 


-1.0(-2) 


-6.1(-2) 


-0.1 


4l4-3o3 


2640.47 


-3.4 


-3.1(-2) 


-0.2 


-0.4 


221-1 10 


2773.97 


-7.5(-3) 


-6.8(-2) 


-0.4 


-0.9 


para-H20 transitions 


2ii-2o2 
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However, due to the lack of more o- and P-H2O observa- 
tions and the poor spectral resolution of the ISO satellite, the 
ISO-LWS data are not decisive on which model is favoured. The 
Herschel mission will bring additional (line profile) informa- 
tion on water transitions with somewhat lower energies, which is 
needed to decide between the different mo dels. Therefore, pr epa- 
ration work, as presented in this paper and lDotv et alj ( l2005h . are 
valuable to prepare for future HIFI observing programs of high- 
mass star-forming regions. 



5.2. Parameter dependency 

5.2.1 . Steepness of density gradient 

We examine the influence of the adopted density and tempera- 
ture distribution, i.e., adopting an or = 2 profile, on the resulting 
line profiles and strengths. In analogy to the models descibed in 
Section |4n we assume the water abundance, ranging from IQ-'*' 
to 10"^, to be uniformly distributed. Resulting line profiles are 
plotted in Fig.|7] integrated brightness temperatures are listed in 
Table |5XT] 

In general, the line shapes are similar to the profiles pre- 
sented in Section \4~T\ However, the 22i-2i2, 3o3-2i2, and the 
ISO transitions are found in absorption, in contrast to the a -I 
models. This is a consequence of a steeper density slope, pro- 
viding more material within the Herschel beam. Note that the 
optical depth in the lines increases by a factor of ~5-10 for 
the high-lying transitions. The line strengths of the ground-state 
lines, i.e., lio-loi and ln-Ooo, are similar to the strengths for the 
a = 1 models. This is not surprising since in the outer regions 
the temperature distribution follows the distribution of the a -I 
model. 

In addition, the effect of an a -2 temperature and density 
profile on the line profiles and strengths is tested for models with 
an outflow cavity, described in Section [43] We find that the line 
shapes are similar to the profiles for the models with an o' = 1 
density gradient. However, the line strengths drop by a factor of 
~ 2 to 3 for the different transitions. 
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Fig. 7: Line profiles of ortho- and para-H20 transitions for mod- 
els with a -2. The water abundance, i.e., 0-H2O + P-H2O, is uni- 
formly distributed across the region and ranges from IQ-'*' to 
10-\ with 10-'" {solid), 10"'' (dotted), 10"^ (dash-dotted), 10"^ 
(dashed). 
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5.2.2. Opening angle 

In this section we briefly discuss the influence of the adopted 
outflow opening angle. We adopt density, temperature and abun- 
dance distributions from Section 14.31 Outflows with a cavity 
half-opening angle of 7.5° and 22.5° are computed. It is found 
that the adopted opening has minor influence on the line shapes. 
However, the line strengths drop (increase) by a factor of ~ 2 - 3 
(2) for models with an half-opening angle of 7.5° (22.5°). 

5.2.3. Outflow water abundance 

In this section we briefly discuss the influence of an increased 
outflow water abundance by a factor of 10, i.e., 10"^. We adopt 
density, temperature and abundance distributions from Section 
14.31 It is found that the adopted outflow water abundance has 
minor influence on the line shapes. However, the line strengths 
increase by a factor of ~ 2 to 10. Note that the lowest gain in 
the line strengths is found for the ground-state transitions, i.e., 
lio-loi, 2i2-l()i, and ln-Ooo, and that the line strengths of the 
higher-lying transitions, e.g., 3i2-3o3, ^n-^ii, 32i-3i2, 221-I01 
and 220-2 n, are increased by a factor up to ~10. This difference 
in enhancement factor can be explained due to the high gas and 
dust temperature in the outflow, favouring the excitation of the 
high-lying transitions. 

6. Summary and Conclusions 

We have constructed models to examine the excitation of wa- 
ter in the high-mass star-forming region AFGL 259 1 . Depending 
on the adopted density, temperature and abundance structure, a 
completely different set of line profiles and strengths is found. 
Hence, the line profiles are very sensitive to the adopted physical 
and chemical structure. We have found that (/) the ground-state 
transitions Iiq-Iqi, 2i2-1oi and ln-Ooo, with relatively low up- 
per energy levels (< 1 10 K), become highly optically thick in the 
outer regions. The line profiles for these transitions, are mainly 
dominated by the emission from the outer regions, and are there- 
fore not useful to put constraints on the water abundance in the 
inner regions. However, (//) the emission from lines with higher 
upper energy levels is dominated by the emission originating in 
the inner regions, and are therefore useful to probe the water 
abundance in the warm inner regions. (///) For models with an 
outflow cavity, the outflow feature (blue peak less strong than 
the red peak) is best seen in the ground-state transitions of o- 
and P-H2O. (iv) The influence of a moderate disk (few 100 AU 
in size) in the centre of the AFGL 2591 region does not change 
the water line profiles and strengths within the Herschel beam. 

The Herschel mission will thus greatly help to understand 
the structure of high-mass protostellar objects, and consequently 
the formation process of high-mass stars. 
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